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Abstract. Nonradial oscillations of relativistic neutron stars with a solid crust are computed in the relativistic Cowling approx- 
imation, in which all metric perturbations are ignored. For the modal analysis, we employ three-component relativistic neutron 
star models with a solid crust, a fluid core, and a fluid ocean. As a measure for the relativistic effects on the oscillation modes, 
we calculate the relative frequency difference defined as 6a ja = (<tgr - a N )/acR, where <tq R and a R are, respectively, the 
relativistic and the Newtonian oscillation frequencies. The relative difference 6a I a takes various values for different oscillation 
modes of the neutron star model, and the value of 6a I a for a given mode depends on the physical properties of the models. 
We find that \6a/a\ is less than ~ 0.1 for most of the oscillation modes we calculate, although there are a few exceptions such 
as the fundamental (nodeless) toroidal torsional modes in the crust, the surface gravity modes confined in the surface ocean, 
and the core gravity modes trapped in the fluid core. We also find that the modal properties, represented by the eigenfunctions, 
are not strongly affected by introducing general relativity. It is however shown that the mode characters of the two interfacial 
modes, associated with the core/crust and crust/ocean interfaces, have been interchanged between the two through an avoided 
crossing when we move from Newtonian dynamics to general relativistic dynamics. 
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1. Introduction 

After the discovery of the r mode instability in neutron stars driven by gravitational radiation reaction (Andersson 1998; Friedman 
& Morsink 1998), nonradial oscillations of relativistic neutron stars has attracted much wider interest than before in astrophysics 
(see a recent review, e.g., by Andersson & Kokkotas 2001). The relativistic formulation of nonradial oscillations of fluid neutron 
stars was first given by Thorne and his collaborators (e.g., Thorne & Campolattaro 1967; Price & Thorne 1969; Thorne 1969a,b), 
and later on extended to the case of neutron stars with a solid crust in their interior (Schumaker & Thorne 1983, see also Finn 
1990). Since one of the main concerns of these studies was gravitational waves generated by the stellar pulsations (e.g., Thorne 
1969a), it was essential to include the perturbations in the metric g a p in order to obtain a consistent description of the gravitational 
waves. However, it was the metric perturbations that made it extremely difficult to treat both analytically and numerically the 
oscillation modes in relativistic stars. 

McDermott, van Horn, & Scholl (1983) introduced a relativistic version of the Cowling approximation, in which all the 
Eulerian metric perturbations 5g a p are neglected in the relativistic oscillation equations derived by Thorne & Campolattaro 
(1967). Under the relativistic Cowling approximation, they calculated the p-, /-, and g- modes of relativistic fluid neutron 
star models to examine how the oscillation modes depend on the model properties. The relativistic Cowling approximation 
employed by McDermott, van Horn, & Scholl (1983) was shown to be good enough to calculate the p-modes for non-rotating 
poly tropic stars by Lindblom & Splinter (1990), who compared the fully relativistic oscillation frequencies to those obtained 
in the relativistic Cowling approximation. Assuming slow rotation, Yoshida & Kojima (1997) have also carried out similar 
computations for the /- and p-modes of polytropic models and they confirmed the good applicability of the relativistic Cowling 
approximation. Quite recently, Yoshida & Lee (2002) have shown in the relativistic Cowling approximation the existence of the 
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relativistic r modes[| with I = m, which are regarded as a counter part of the Newtonian r modes. We believe that the relativistic 
Cowling approximation is quite useful to investigate the oscillation modes of relativistic stars, although we understand that under 
the approximation we cannot discuss inherently relativistic oscillation modes like w-modes (see, e.g., Andersson, Kojima, & 
Kokkotas 1996). 

In their Newtonian calculations, McDermott, van Horn, & Hansen (1988) have shown that cold neutron stars with a solid 
crust can support a rich variety of oscillation modes (see also Strohmayer 1991; Lee & Strohmayer 1996; Yoshida & Lee 2001). 
The purposes of this paper are to calculate the various oscillation modes of relativistic neutron stars with a solid crust in the 
relativistic Cowling approximation, and to discuss the effects of general relativity on the modal property of the oscillation modes. 
We regard this paper as an extension of the studies by McDermott, van Horn, & Scholl (1983) and McDermott, van Horn, & 
Hansen (1988). In §2, we derive relativistic oscillation equations for the solid crust in the relativistic Cowling approximation 
following the formulation developed by Schumaker & Thorne (1983) and Finn (1990). Numerical results are discussed in §3 for 
nonradial modes of three relativistic neutron star models with a solid crust. §4 is devoted to discussion and conclusion. In this 
paper, we use units in which c - G - 1, where c and G denote the velocity of light and the gravitational constant, respectively. 



2. Formulation 

For modal analysis of neutron stars with a solid crust, we solve general relativistic pulsation equations derived under the rela- 
tivistic Cowling approximation, in which all the metric perturbations in the matter equations are neglected. We assume that the 
solid crust is in a strain-free state in the equilibrium unperturbed state and the strain in the crust are generated by small amplitude 
perturbations superposed on the unperturbed state (see, e.g., Aki & Richards 1980). For the background unperturbed state, it is 
therefore possible to assume a static and spherical symmetric state, for which the geometry is given by the following line element: 

(ds 2 \ 0) = g ( o )a /idx a dx (i = -e 2v(r) dt 2 + e 2A{r) dr 2 + r 2 d6 2 + r 2 sin 2 6dip 2 , (1) 

where the subscript (0) is used to refer to the unperturbed state. The four-velocity of the stellar material in the equilibrium is 
given by 

u? 0) = e-^f, (2) 

where f denotes the timelike Killing vector of the unperturbed spacetime. 

In the relativistic Cowling approximation, the basic equations for pulsations are obtained from the energy and momentum 
conservation laws: 

u a VpTa = , (energy conservation law) (3) 
q" VgTf = , (momentum conservation law) (4) 

where V„ is the covariant derivative associated with the metric, is the energy-momentum tensor, and is the projection 
tensor with respect to the fluid four-velocity u a , which is defined by 

$ = (5) 

where a a denotes the Kronecker delta. 

For relativistic analysis of vibration of the solid crust in neutron stars, we follow the formulation given by Schumaker & 
Thorne (1983) and Finn (1990). If we define the rate of shear (j a p and the rate of expansion 6 as 

o~a/3 = ^ qlq S p (V 7 m 5 + V s u y ) - ^ , (6) 

= (f a Vpu a = V a if , (7) 

the shear strain tensor Y, a p is determined, in terms of <r a p and 6, as a solution of the differential equation given by 
_ 2 

L u ^aj3 — g #£oy8 + °~a/3 > (8) 

where L u is the Lie derivative along the four-velocity field of the matter (Carter & Quintana 1972). Once equation (|j|) is solved to 
give the shear strain tensor "Z a p in terms of <r a p and 6 (see Schumaker & Thorne 1983, Finn 1990), the total stress-energy tensor 
is given by 

T a p = p U a U/3 + pq a p- 2/i I, al 3 , (9) 

1 We should note that lots of studies about the general relativistic r-mode have been done so far (e.g., Kojima 1998; Kojima & Hosonuma 
2000; Lockitch, Andersson & Friedman 2001; Yoshida 2001; Ruoff & Kokkotas 2001). 



S. Yoshida & U. Lee: Nonradial Oscillations of Neutron Stars 



3 



where p stands for the isotropic shear modulus, and we have assumed a Hookean relationship between the shear strain and stress 
tensors. Here, p and p mean the mass-energy density and the isotropic pressure, respectively. Note that T.(0) a /3 = since we have 
assumed the strain-free state in the equilibrium of the star. 

To discuss nonradial oscillations of a star, we employ a Lagrangian perturbation formalism (see, e.g., Friedman & Schutz 
1975; Friedman 1978), in which a Lagrangian displacement vector is introduced to connect fluid elements in the equilibrium 
state to the corresponding elements in the perturbed state. In this formalism, the Lagrangian change AQ in a quantity is related 
to its Eulerian change 5Q by 

AQ = 6Q + L C Q, (10) 

where denotes the Lie derivative along the displacement vector g*. If we apply this formalism to oscillations of relativistic 
stars, we have, for example, 

Agafi = Sgap + + V^ a , (11) 

and 

Au a = X -u a \fu y Agfsy. (12) 
The Eulerian perturbation in the velocity field 5$ = q^dif is then given by 

8u a = q m a p {L u ^f . (13) 

Notice that 8g a p = and hence 6u a = 8u a in the relativistic Cowling approximation. Because the background state is in a 
hydrostatic equilibrium, the time dependence of all the perturbed quantities can be given by e" 7 ', where cr is a constant frequency 
measured by an inertial observer at the spatial infinity. The relation between the Lagrangian displacement C, a and the velocity 
perturbation 8u" is then given as an algebraic equation: 

Su a = io-e- vir) C- (14) 

Note that the gauge freedom in ( a has been used to demand the relation u a ( a = 0. The velocity field of a neutron star in a 
perturbed state may be written as 

u a = u? ) + rj8u a + 0(tj 2 ), (15) 

where 77 is a small expansion parameter introduced for our convenience. The mass-energy density p and the pressure p in the 
perturbed state are also given by 

lYSp \ n 

P = P(0) + n (p ( o) + P(p)) C A a + 0(rf) , (16) 

\P(P) I 

p= P{Q) +i 1 8p + 0{T 1 2 ), (17) 
where A a is the relativistic Schwarzschild discriminant defined as 

A a = V a p (0) - -J— V a p( ) , (18) 

P(0) + P(0) 1 P(Q) 

and F is the adiabatic index defined as 

r _ P(Q) +P(Q) ( dP(Q) \ (19) 
To derive equation (16), we have used the adiabatic condition for the perturbations: 

Ap = Tm Ap. (20) 

P(0) + P(0) 

Assuming spherical symmetry of the background spacetime, we may expand the perturbations in terms of appropriate tensor 
spherical harmonic functions. The Lagrangian displacement, £ k and the pressure perturbation, 8p/(p(0) + p(0)) can be written as 

C = rSi(r)Y?<<9,<p)e UTt , (21) 

. I dY™(6,<p) 1 dY?{6, <p)\ . , 

\ o6 sin0 dip j 
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1 / dY"'(e,tp) dY?{6,(p) 
— \H[(r) — — — + 7>(/")sin0 — '— Y 
sin 2 6» \ dip 80 



I* = —j- ( H,(r) l J + 7>(r) sin 6 r J ) e im , (23) 



<V ' = 5Ui{r)Y?{6,<p)e i<Tt , (24) 



P(0) + P(0) 



where Yf is a spherical harmonic function (Regge & Wheeler 1957; Thorne 1980). 

By substituting equations ([)]), (15), (16), and (17), together with the perturbations defined by equations (21) to (24), into 
equations (||) and (Q), and collecting the terms proportional to the small expansion parameter 77, we obtain a system of perturbation 
equations for the solid crust: 

dr \ a-i J 

= \(-3-U2 + Ui-e 2A c 1 d- 2 )Vi + —(3a2 + 2a 1 )\zi+tv 2 -4—)z2 
dr { a 3 ) \ a-i) 



2c?2 



+ |Vi-2ai + f(/+l)*S +«"/(/+ l)z«, (26) 

dz% 1 3 e 2A 

r^ = -e 2A Zl + —Z4, (27) 
dr a\ 

r — = - -Vi + 6r — zi z 2 

dr \ QT3 / 

- |c 1( r 2 yi + 2ori - (05 + or 3 ) /(/ + 1)J Z3 - (3 + U 2 - V 2 )Z4 , (28) 
dzs e 2A 

r-^ = —ze, (29) 
dr a\ 



dZ6 

r 

dr 



(3 + U 2 - V 2 )z 6 - (c^V, - <*!(/' + 1)(/' - 2)} z 5 , (30) 
where the dependent variables z\ to Z6 are defined as 

zi=St(r), (31) 

z 2 = 2axe- A j r (re A S,(r)) + |r - | a, J |^ i. (rVS,(r)) - /(/ + l)#/(r)J , (32) 

Z3=H t (r), (33) 

Z4 = ffl ( e -V^ +5/W ), (34) 

Z5 = r?(r), (35) 
-zj dT v (r) 

z 6 = a x e r — - — (36) 
dr 

and the various quantities which appear in the coefficients are 

ai = — , a 2 =r-|ari, ar 3 =r+^ori, (37) 
P(0) 3 3 

+ V 2 = ^r^, C/^f^r^fr^, U 2 = r d 4, (38) 



££>(0) / rfr S P(0) dr \drj dr\ drj ' dr 
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ei = 8 *re-*>m , (39) 



R 3 \dr 

eU = / 1 _ g ^£W\ _ A/(/ . )= ■ 4rTrpi(h(ll ^ (40.) 



Jo 



^ ?i /„ ? M(r)\ dA 11 [ * 7 M(r)\ 

r - = ee 2 * + -iij , r - = ^ (Wp (0) - -iij . (41) 

Here, M = M(/?) and /? are the mass and the radius of the star, and a- - cr -\/R 3 /M is the dimensionless frequency. In equations 
(|3~8|)-([4l]), we have introduced a factor e to indicate the hidden factor 1/c 2 that represents the strength of the general relativistic 
degree. Fully relativistic oscillation equations for the solid crust are regained for s = 1, and the Newtonian oscillation equations 
are obtained in the limit of s — > 0, which are the same as those derived by McDermott, van Horn, & Hansen (1988). We can see 



that equations (|25|)-(|28|) are decoupled from equations (29)-(p0|) because the background spacetime is spherical symmetric. The 
former describe spheroidal (or polar parity) oscillations, and the latter toroidal (or axial parity) ones. 
In the fluid regions, the oscillation equations to be solved are given by 



dyi 
dr 



r d 2L = (e 2A Cl d- 2 + rA r )y x - (U + rA r )y 2 , (43) 
dr 

where 

d Y 1 

— I SUM = c^H,. (44) 

Note that equations (42) and (43), which come from equations (^5|)-(|28|) in the limit of jj. — » 0, are essentially the same as those 
derived by McDermott, van Horn, & Scholl (1983) for spheroidal oscillation modes. Equations (|29"|)-(|3"o|), which describe toroidal 
torsional oscillations, become trivial in the limit of ft — > because there is no restoring force that is responsible for the toroidal 
torsional oscillations. To calculate the spheroidal modes, equations (^5|)-(^8|) are integrated in the solid crust, and equations (^2|)- 
( 43 ) in the fluid core and the surface ocean. For the toroidal torsional modes, equations (^9|)-(^o|) should be solved only in the 



solid crust. 

For the spheroidal modes the outer boundary condition is given at the stellar surface (r = R) by Ap = 0, which reduces to 

yi-y 2 = 0, (45) 

and the inner boundary condition is the regularity condition at the stellar center given by 

c 1 (7 2 y 1 -ly 2 = Q. (46) 

The jump conditions at the interfaces between the fluid and the solid regions are given by the continuity conditions, across the 
perturbed interface, of the stress P a on the interface, where 

P a = {pq a p-2 l i\Z ap )N p , (47) 

and is the unit vector normal to the perturbed interface (see, e.g., Finn 1990). The normal form of the core/crust or crust/ocean 
interface is given by d(r — r]^ r ) = riadx" for the radial coordinate of the perturbed interface given by r — Rj + rj £ r , where is the 
radius of the core/crust or crust/ocean interface in the unperturbed state. Thus, the corresponding unit normal one-form have 
the components as follows: 

N, = -e A T 1 io-rSi(r)Y';X0,ip)e i ' T ' + 0(tj 2 ) , 

N r = e A + 0(i] 2 ) , (48) 
N e = -e A jj rS ,(r) l — - e^' + 0( V 2 ) , 

O0 

N v = -e A rS,(r) 1 Y e m + 0(t, 2 ) , 

The jump conditions for the spheroidal modes at the interface are then given by 



yi = zi , ViCyi - yi> - zi , za = 



(49) 
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For the toroidal torsional modes, the boundary conditions are applied at the top and the bottom of the solid crust, and they are 

ze = 0. (50) 

By imposing the boundary conditions given above, we can solve our basic equations as an eigenvalue problem with respect to 
the eigenvalue &. Here, we employ a Henyey type relaxation method to obtain numerical solutions to our basic equations (see, 
e.g., Unno et al. 1989). 



3. Numerical Results 

Neutron star models that we use in this paper are the same as those used in the modal analysis by McDermott, Van Horn, & 
Hansen (1988). These models are taken from the evolutionary sequences for cooling neutron stars calculated by Richardson et al. 
(1982), where the envelope structure is constructed by following Gudmundsson, Pethick & Epstein (1983). They are composed 
of a fluid core, a solid crust and a surface fluid ocean. The interior temperature is finite and is not constant as a function of the 
radial distance r. The models are not barotropic and the Schwarzschild discriminant \A\ has finite values in the interior of the 
star. The models we use are called NS05T7, NS05T8, and NS13T8, and their physical properties such as the total mass M, the 
radius R, the central density p c , the central temperature T c and the relativistic factor GM/c 2 R are summarized in Table 1 (for 
other quantities, see McDermott et al. 1988). 

To classify the various oscillation modes of the three component models, we use almost the same nomenclature as that 
employed by McDermott et al. (1988). We let pk refer to the acoustic modes of the k th overtone. The internal gravity modes 
confined in the surface ocean are denoted as g s k and those in the fluid core as g c k , where k indicates the overtone number. The 
eigenfrequencies of the / modes are usually found between the p\ and g\ modes. Associated with the fluid-solid interfaces in the 
models, there are two interfacial modes, which we denote as i\ and z'2 such that cr{i\) < cr(z 2 ). The Sk modes are spheroidal shear 
dominated modes of the k th overtone, the amplitudes of which are strongly confined in the solid crust. The t k modes are toroidal 
shear dominated modes of the k th overtone propagating only in the solid crust. Note that the p^, gk-, /-, z'i( 2 ), and Sk modes are 
classified as spheroidal modes while the tk modes as toroidal modes. 

In the Newtonian limit of s — > 0, we calculate various oscillation modes with I = 2 for the three neutron star models NS05T7, 
NS05T8, and NS13T8, and tabulate the Newtonian eigenfrequencies, which we denote as cr N , in Tables 2 through 4. These 
tables confirm that the Newtonian frequencies obtained in this paper are in good agreement with those computed by McDermott 
et al. (1988). Assuming s = 1, we also compute the corresponding relativistic oscillation modes with I = 2, and tabulate the 
relativistic eigenfrequencies, which we denote as cr CR , m Tables 2 through 4, in which the relative frequency differences defined 
as bcrjcr = (cr CR - cr a?)/ cr q R are also given. The relative difference Scr/cr takes various values for different oscillation modes, and 
the value of dcr/cr of a given mode is also dependent on the physical properties of the neutron star models. We find \6cr / cr\ <0.1 
for most of the oscillation modes we calculate. The elastic Sk and tk modes with k > 1 have particularly small values of \6cr/cr\, 
which are at most a few percent. However, there are some exceptions, and the values of 6cr/cr become as large as -0.45 for 
the to modes and ~ — 1 for the g s - and g c -modes for the model NS13T8, which is the most compact one among the three. If 
we consider the gravitational redshift as one of the general relativistic effects that are responsible for the difference from the 
Newtonian oscillation frequency, we may have &cr = e v(f) cr N and hence (cr CR - cr N )/cf CR = 1 - e~ v{r \ which is negative since 
v(r) < 0. Although this estimation gives a consistent result for the surface g s modes confined in the surface ocean and for the core 
g c modes trapped in the fluid core, it does not always give good estimations of the difference as suggested by the tables. In fact, 
we have positive dcrjcr for some oscillation modes. This is understandable, however, since the eigenfunctions usually possess 
large radial extent in the interior of stars, and the eigenfrequencies are dependent on the equilibrium quantities which contain 
terms and factors which are absent in the Newtonian oscillation equations. 

For the various oscillation modes of NS13T8, we display the displacement vector £ a versus the radial coordinates r/R or 
log(l - r/R) in Figures 1 through 8, where the amplitude normalization of the eigenfunctions is given by y\ = 1 at r — R for the 
spheroidal modes and zs — 1 at r — R out for the toroidal modes with R out being the radius of the ocean-crust interface. Notice 
that for the g c modes, the normalization condition //;(/?,„) = 1 is adapted because the modes are very insensitive to the behavior 
of the eigenfunctions in the surface ocean, where is the radius at the bottom of the crust. The eigenfunctions obtained in 
the Newtonian calculations are also plotted for comparison in each of the figures. From the figures 1 through 8, we can see that 
the basic properties of the eigenfunctions of the relativistic modes are not very much different from those of the corresponding 
Newtonian modes, except for the interfacial modes (see the next paragraph). In this sense, we can say that general relativity does 
not bring about any essential changes in the modal properties represented by the eigenfunctions. 

As shown by Figures 3 and 4, the relativistic eigenfunctions of the interfacial modes i\ and z'2 are at first sight quite different 
from the Newtonian eigenfunctions. Since the eigenfunctions of the relativistic z'i(2) mode are rather similar to those of the 
Newtonian ?2(i) mode, it is tempting to make a guess that the modal properties have been exchanged between the two modes 
when we move from Newtonian dynamics to relativistic dynamics. In their Newtonian calculation, McDermott et al. (1988) 
showed that the properties of the z modes are very sensitive to the changes in the physical quantities near the interfaces, and 
mentioned a numerical experiment in which they computed the two interfacial modes by artificially reducing the bulk modulus 
/j. and found an avoided crossing between the two modes through which the modal characters are interchanged with each other. 
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Inspired by this report, we have carried out a numerical experiment, calculating the two interfacial modes as a function of s for 
the model NS13T8. Plotting the frequencies of the two modes versus s in Figure 9, we confirm the occurrence of an avoided 
crossing between the two modes at s ~ 0.67, through which the modal characters of the two have been exchanged with each 
other. As s increases from s - (Newtonian limit), the effective local shear modulus fi in the solid crust decreases because of the 
gravitational redshift effects, as a result of which the two interfacial modes experience the avoided crossing when we move from 
Newtonian dynamics to relativistic dynamics. 

4. Conclusion 

In this paper, we have calculated a variety of oscillation modes of relativistic neutron stars with a solid crust in the relativistic 
Cowling approximation, in which all metric perturbations are ignored. We find \6cr/o-\ is less than ~ 0. 1 for most of the oscillation 
modes with I = 2, although there are some exceptions such as the surface gravity modes, the core gravity modes, and the nodeless 
toroidal torsional modes. We also find that the essential modal properties represented by the eigenfunctions are not strongly 
affected by introducing general relativity in the sense that the relativistic eigenfunctions of a mode have almost the same properties 
as the corresponding Newtonian eigenfunctions. An exception may be the two interfacial modes whose eigenfunctions differ from 
their Newtonian eigenfunctions, and we have shown that this can be explained as a result of the mode exchange between the two 
modes through an avoided crossing that occurs when we move from Newtonian dynamics to relativistic dynamics. 

Needless to say, nonradial oscillations of neutron stars should be treated within the framework of general relativity because 
of their strong gravity. At present, however, we cannot investigate the pulsations of neutron stars in all the aspects by using fully 
general relativistic formalism because of the complexity general relativity brings about.^j One of the simplest ways to study the 
oscillations of a neutron star is to treat the problems within the framework of Newtonian dynamics, that is, to solve the Newtonian 
oscillation equations for neutron star models constructed with the Newtonian hydrostatic equations. However, it is obvious that 
this Newtonian treatment cannot be fully approved for neutron stars since the Newtonian equilibrium models do not give us 
correct mass and radius for the stars. To reduce the distance between fully relativistic and Newtonian calculations of nonradial 
oscillations of the stars, we may use relativistic equilibrium models and solve the Newtonian oscillation equations, which was the 
strategy taken by McDermott, van Horn, & Hansen (1988). In this paper, as an extension of the study by McDermott, van Horn, 
& Hansen (1988), we solved relativistic oscillation equations derived in the relativistic Cowling approximation for relativistic 
neutron star models with a solid crust. Our calculation confirmed that the classification scheme McDermott, van Horn, & Hansen 
(1988) employed for the oscillation modes of neutron stars with a solid crust is valid even if we integrate relativistic oscillation 
equations. 

Acknowledgements. S.Y. acknowledges financial support from the Portuguese FCT through a Sapiens project, number 36280/99. 
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Table 1. Neutron Star Models 



Model 


M(M ) 


R(km) 


Pc (g cm 3 ) 


r c (K) 


GMKc'R) 


NS05T7 


0.503 


9.839 


9.44 x 10 14 


1.03 x 10 7 


7.54 x 10~ 2 


NS05T8 


0.503 


9.785 


9.44 x 10 14 


9.76 x 10 7 


7.59 x 10~ 2 


NS13T8 


1.326 


7.853 


3.63 x 10 15 


1.05 x 10 8 


2.49 x 10-' 



Table 2. Eigenfrequencies & (I = 2) of NS05T7 



mode 




(TOR 


Scr/cr 


(S pheroidal) 








s\ 


1.998 x 10" 5 


1.625 x 10- 5 


-0.230 


Si 


2.037 x 10~ 3 


1.760 x 10~ 3 


-0.158 




3.230 x 10~ 3 


2.787 x 10- 3 


-0.159 


h 


8.282 x 10- 3 


7.499 x 10" 3 


-0.104 


h 


1.029 x 10"' 


9.717 x 10" 2 


-0.059 




3.093 x 10-' 


3.159 x 10-' 


0.021 




5.556 x 10-' 


5.548 x 10-' 


-0.001 


f 


1.886 x 10° 


1.823 x 10° 


-0.034 


Pi 


4.038 x 10° 


3.978 x 10° 


-0.015 


P2 


4.736 x 10° 


4.710 x 10° 


-0.006 


(Toroidal) 








to 


4.019 x 10- 2 


3.619 x 10- 2 


-0.111 


h 


3.286 x 10-' 


3.279 x 10~' 


-0.002 


h 


5.633 x 10-' 


5.609 x 10-' 


-0.004 


h 


7.346 x 10- 1 


7.355 x 10-' 


0.001 



Table 3. Eigenfrequencies <x (7 = 2) of NS05T8 



mode 




0~GR 


Scr/cr 


(5 pheroidal) 








s\ 


2.248 x 10- 4 


1.821 x 10- 4 


-0.234 


Si 


1.297 x 10~ 2 


1.102 x 10- 2 


-0.177 


8\ 


1.507 x 10- 2 


1.289 x 10- 2 


-0.169 


h 


5.346 x 10- 2 


5.894 x 10- 2 


0.093 


h 


1.002 x 10"' 


9.236 x 10- 2 


-0.085 




3.251 x 10-' 


3.299 x 10-' 


0.014 


S2 


5.755 x 10-' 


5.746 x 10-' 


-0.002 


f 


1.873 x 10° 


1.811 x 10° 


-0.035 


Pi 


4.102 x 10° 


4.027 x 10° 


-0.019 


P2 


4.766 x 10° 


4.747 x 10° 


-0.004 


(Toroidal) 








to 


3.999 x 10~ 2 


3.600 x 10~ 2 


-0.111 


h 


3.419 x 10-' 


3.410 x 10-' 


-0.003 


t2 


5.841 x 10-' 


5.812 x 10-' 


-0.005 


h 


8.148 x 10~' 


8.128 x 10-' 


-0.002 
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Table 4. Eigenfrequencies <x (/ = 2) of NS13T8 



mode 


0~N 


0~GR 


So~jo~ 


(S pheroidal) 








s c 

a\ 


1.269 x 10~ 4 


6.017 x 10~ 5 


-1.109 


62 


6.115 x 10~ 3 


3.126 x 10~ 3 


-0.956 


2 s 
6] 


7.548 x 10~ 3 


3.905 x 10~ 3 


-0.933 


h 


1.862 x 10~ 2 


2.357 x 10~ 2 


0.210 


h 


3.326 x 10~ 2 


2.985 x 10~ 2 


-0.114 


Si 


4.512 x 10~' 


4.528 x 10~' 


0.004 


s 2 


7.704 x 10~' 


7.685 x 10~' 


-0.002 


f 

J 


1.434 x 10° 


1.249 x 10° 


-0.148 


Pi 


3.967 x 10° 


3.038 x 10° 


-0.306 


P2 


5.511 x 10° 


4.602 x 10° 


-0.197 


(Toroidal) 








to 


1.897 x 10~ 2 


1.309 x 10- 2 


-0.449 


h 


4.534 x 10-' 


4.533 x 10-' 


-0.000 


h 


7.727 x 10-' 


7.695 x 10-' 


-0.004 


h 


1.105 x 10° 


1.100 x 10° 


-0.004 



■S 0.5 
3 



s 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — IT 



r S/R 

r S/R (Newton) 

r H/R 

_.r H/R (Newton) 




J I I I I I I I I I I I I I I I I I L 



0.2 0.4 0. 

r/R 



0.8 1 



Fig. 1. Displacement vectors rSi/R and rfy/R of the g c x mode for the model NS13T8, given as a function of r/R. Here, the 
normalization of the eigenfunction is chosen as H2(Rm) = 1, where /?,„ is the radius at the bottom of the crust. The Newtonian 
eigenfunctions are shown as well as the relativistic ones. 
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„ r S/R 

6 ~ \ r S/R (Newton) 

\ r H/R xlO- 4 

XJ 4 - \ _ r H/R xlO-" (Newtonf 



0, 



E 2 - 




_ g I i i i I i i i I i i i I i i i 

-8 -6 -4 -2 

\og{l-r/R) 

Fig. 2. Displacement vectors rS2/R and rHilR of the g\ mode for the model NS13T8, given as a function of log(l - r/R). Here, 
the normalization of the eigenfunction is chosen as S ziR) = 1 • The Newtonian eigenfunctions are shown as well as the relativistic 
ones. 



30 i i i g — i — i — i — | — i — i — i — | — i — i — r 

'N 

\ 

\ 

20 - \ 



0) 

-a 10 

E o \=-- 



10 - 



I 

i y 



-r S/R 
r S/R xO.l (Newton 
r Hg/R xlO' 2 
r H/R xlO' 2 (Newton) 




-20 



_l I I I L 



-6 -4 -2 

log(J-rA) 



Fig. 3. Same as Figure 2 but for the i\ mode. 
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CD 

Q, 

E 

^-10 
-15 

-20 




r S/R 

r S/R (Newton) 

r Hg/R xW-* 

r H/R xW-? (Newton) 



-6 -4 
log(l -r/R) 



Fig. 5. Same as Figure 4 but for the s\ mode. 
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Fig. 6. Displacement vectors rS2/R and rH-ilR of the / mode for the model NS13T8, given as a function of r/R. Here, the 
normalization of the eigenfunction is chosen as S2W = 1- The Newtonian eigenfunctions are shown as well as the relativistic 
ones. 
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1 

0.8 
0.6 
0.4 
0.2 


-0.2 

0.94 0.96 0.98 1 

r/R 

Fig. 8. Displacement vector rTz/R of the to and the fi modes of the neutron star model NS13T8, given as a function of r/R. Here, 
normalization of the eigenfunction is chosen as T 2 (R„ui) = 1, where R out is the radius at the crust/ocean interface. The Newtonian 
eigenfunctions are shown as well as the relativistic ones. Note that two curves for rTtjR nearly overlap each other. 



,/ 

T L 



_i i i_ 



t (Newton) 
t, 

t t (Newton) 
j i i_ 



i i i | i i i ^ i i i | i i i j i i r 




0.2 0.4 0.6 0.8 1 

£ 



Fig. 9. The avoided crossing between the i\ and i 2 mode of the model NS13T8, where s is the parameter that represents the 
strength of general relativistic effects. The avoided crossing happens at s ~ 0.67. 



